Domain walls in chiral p-wave superconductors: Quasiparticle spectrum and dynamics 
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We calculate microscopically the viscous friction coefficient and the effective mass of domain walls 
separating regions of opposite chirality in p-wave superconductors with ± iky order parameter. 
The domain wall viscosity and inertia are determined by the transitions between different Bogoliubov 
quasiparticle states induced by the domain wall motion. As a by-product, we present a detailed 
analysis of the quasiparticle spectrum, both bound and scattering, in the presence of a general 
^— ■ domain wall with an arbitrary phase difference between the domains. 

(N . 

^ ' PACS numbers: 74.20.Rp, 74.20.-z 
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The properties of topological superconductors and superfluids have attracted a considerable interest recently. The 
defining feature of these systems is that, while the fermionic excitations in the bulk are fully gapped, nontrivial 
^ ' topology of the order parameter manifests itself in the existence of gapless quasiparticles localized near the order 
parameter inhomogeneities, such as Abrikosov vortices, domain walls, or the sample boundaries. One of the most 
r ' studied examples is the chiral p-wave spin-triplet state, whose experimental realizations include the superconducting 
state of Sr2Ru04 (Ref. or thin films of superfluid '^He-A (Ref. '3). The chiral p-wave state in fermionic superfluids 
^ ' is closely related to the Moore-Read Pfaffian quantum Hall state. ^ 

c/5 , The momentum-space order parameter of a chiral p-wave superconductor is proportional to ± iky. In the absence 
of external magnetic field, it is two- fold degenerate: the states k^ + iky and kx — iky, which are obtained from each 
other by time reversal, have the same energy. Therefore, superconducting states with opposite chiralities separated 
by domain walls (DW) might form in different parts of the system. Indeed, there is evidence of the DW existence in 
Sr2Ru04 (Refs. iandi) and also in slabs of superfluid ^He (Ref. @)- The DW formation costs gradient energy due 
^ I to the spatial variation of the order parameter. In contrast to ferromagnets, which break up into domains in order 
Q ' to minimize the net magnetic moment, there is no similarly compelling energy reason in a neutral superfluid. One 
O I possible mechanism is that domains are spontaneously formed upon cooling across the phase transition due to the 
sample inhomogeneity. Alternatively, an increase in the gradient energy might be compensated by the creation of 
I ' low-energy quasiparticles bound to the DW, which is particularly effective in one-dimensional systems.™ 
• While the static properties of the DWs in various realizations of the chiral p-wave state have been extensively 
Qf^ ' studied^'"— their dynamics has received comparatively little theoretical attention, see Ref. [l^. The motion of a 
0^ • different type of planar defects, namely, an interface between the A and B phases of superfluid ^He was studied in 
Refs. [10 and [U, where it was pointed out that the scattering of Bogoliubov quasiparticles by the moving interface 
results in an effective friction force. Similar ideas can be applied in our case as well. The DW, which is assumed 
to be moving uniformly as a whole, suffers viscous friction and acquires mass due to its interaction with fermionic 
quasiparticles, see Sec. [ill To obtain the DW dynamic characteristics we employ the effective bosonic action formalism, 
see Sec. IIIIl The Gaussian effective action for the DW essentially depends on the Bogoliubov quasiparticle spectrum, 
both bound and scattering states, in the presence of a static DW. The latter is studied in detail for a general DW 
structure in Sec. IIVI Although some bits and pieces about the properties of the DW quasiparticle spectrum can 
be found scattered in the literature, we believe it is useful to present a complete picture in one place. Some of the 
more technical details are discussed in four appendices. Finally, the DW friction coefficient and the zero-temperature 



' effective mass are calculated in Sec. |Vl Throughout the paper we use the units in which h = ks = c = 1. 



II. THE MODEL 



Let us consider a two-dimensional triplet p-wave fermionic superfluid or superconductor. The gap function is a spin 
matrix given by i{&d)a2, where d = z{r]ikx+rj2ky) / kp describes triplet pairing, kp is the Fermi wavevector, and <t are 
the Pauli matricesj22, The order parameter r] = (771, 772) is characterized by two planar components, which can depend 
on coordinates and, in the dynamic case, time. We assume that the band dispersion is isotropic: ^(fc) = {k'^ — kp) / 2m* , 
with the effective mass m* (our results can be straightforwardly generalized for the case of anisotropic dispersion). 
We also neglect disorder, as well as the effects related to the electric charges, such as the Meissner screening of the 
external or internal magnetic fields. 
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For a static DW, one can choose the x-axis along the normal and write the order parameter as 

»? = (hi 1,1^2 |e-^>'^ (1) 

where the amplitudes of the components, the relative phase 7, and the common phase (p all depend on x. Assuming that 
the most stable superconducting state in the bulk is described by one of the two degenerate chiral states rj oc (1, ±i) 
and allowing for a nonzero phase difference between the two domains, we have the following expression for the order 
parameter asymptotics far from the DW: 

77(0;) = Ao(l,i), 

J7(a;) = Aoe*x(l,-i): .t ^ +00. ^ > 

Here x is a parameter which depends on the microscopic details (0 < X < tt). Its value is fixed by the condition of 
vanishing supercurrent across the DW, see a discussion of this point in Appendix [A] 

An exact analytical expression for the DW structure is not available and a variety of approximations have been 
proposed in the literature. For instance, the amplitudes of both components can be put constant: [771 1 = [772 1 — Aq 
(Ref. Q)- Alternatively, one can assume constant phases: r] — Ao(l, i/) or rj = Ao(«/, 1), where a real function 
f{x) varies between 1 at a: = —00 and —1 at a; = 00 (Refs. l^lsl llOl. and fill). Other possibilities include T7 — 
Ao(cos 6, i sin9), where Q{x) varies between and tt (Refs. [H and IITI). and (?7+,7?-) = Ao(e"'^+ cosCiC**^- sin^), 
where ri± — (771 =F 2772 )/v^, 4>± are parameters, and C,{x) varies between and 7r/2 (Refs. [13 and ITsl). In all cases, 
the DW order parameter variation occurs within a region of width around the origin x = 0, with being the DW 
thickness. The precise way in which the DW order parameter varies between the asymptotics given by Eq. ^ is not 
important for our purposes. 

We are interested in the motion of the DW as a whole. Such a motion can be caused, for instance, by the (extremely 
weak) interaction of an external magnetic field with the orbital moment of the Cooper pairsi^ The direction of the 
latter is given by the unit vector I = i{r]* x T])/\'r]* x t]\, which takes opposite values in the two domains [according 
to Eq. ([T]), Z = z sign (sin 7)], thus creating a transverse force on the DW. For a small driving force one can expect 
a linear relation between the DW velocity and the force. According to the fundamental principles of nonequilibrium 
statistical mechanics, one can express linear-response kinetic coefficients, in particular, the DW viscous friction, in 
terms of the equilibrium fluctuation properties. It is legitimate, therefore, to use the Matsubara formalism with the 
magnetic field and the driving force both set to zero. In the absence of external field, the controversy about the 
magnitude of the orbital moment and the lack of its local definition, see Refs. 22^2! and [25l. is not relevant for our 
problem. 

In the bosonic effective action formalism,— the order parameter becomes a two-component dynamical field 'r]{r, r), 
which depends not only on coordinates, but also on the Matsubara time r. We use the following ansatz for a moving 
DW: 

rj{r,T)^ri[x-X{T)], (3) 

where ri{x) is the mean- field order parameter of a static DW and X{t) is the macroscopic coordinate describing a 
uniform displacement of the DW. The effective action can be expressed in terms of X as follows: 

Seff[X] = 5o +T^/C(7/„)X(z/„)X(-z.™), (4) 

m 

where is the mean-field action for the static DW, the second term is the DW dynamic action in the Gaussian 
approximation, and Vm = ^mirT is the bosonic Matsubara frequency. Due to the translational invariance of the 
system, a static, i.e. time-independent, displacement would not affect the action, therefore, ]C{i'm = 0) = 0. 

One can view the DW as a macroscopic quantum object interacting with the equilibrium bath of fermionic quasi- 
particles. Its motion induces transitions between the quasiparticle states, which result in the effective DW friction 
and inertia. If the leading term in the frequency expansion of the kernel /C(j^,„) is given by »7|fm|j then rj can be 
interpreted, in the spirit of the Caldeira-Leggett model, as the viscous friction coefficient, while a term quadratic 
in Vm yields the inertial mass of the DW. The frequency is assumed to be small compared to the gap amplitude, so 
that the DW order parameter profile is not deformed during its motion. Note that there are different ways to define 
the effective mass of a topological defect in superconductor and superfluids, discussed mostly in the context of the 
Abrikosov vortex dynamics, see Ref. H^. Our approach, based on the Matsubara effective action, is similar to the one 
developed for the vortex dynamics in Ref. '29. In general, the dynamics of any stable inhomogeneous order parameter 
structure can be analyzed in this way.— 
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III. DERIVATION OF THE EFFECTIVE ACTION 



Our investigation of the DW dynamics is based on the bosonic effective action for a chiral p-wave superconductor. 
The standard procedure, which involves integrating out the fermionic degrees of freedom^ yields the following 
expression for the action: 

^'"^vJ!'^^ J dr\vir,r)\^-lTT\ng-\ (5) 
where 1/ > is the coupling constant of the triplet p-wave pairing channel, /3 is the inverse temperature, and 

t;r^ 'T'^^)- (6) 



-At (r,T) 

Here ^ = ^(fc), A(r,T) — ai[rji{r,T)kx + ri2{r,T)ky]/kF, and fc = — iV. For a DW moving as a whole, see Eq. ([3|), 
we have A(r, r) — Ao[x — X{t)], where 

A / ^ - mix)k +mix)ky 

Ao(x)=CTi (7) 

kp 

corresponds to the static DW, with the order parameter components given by Eq. ([1]). 

The first term in the effective action ([5]) does not depend on X, while the second one can be expanded in powers 
of the displacement, using ^ Q^^ — E, where 

^o(n,r2;a;„) = 2^-— — — (8) 

is the Green's function at X — 0. Here w„ = (2rt + l)7rT is the fermionic Matsubara frequency, \j) and Ej are the 
eigenfunctions and eigenvalues of the 4x4 Bogoliubov-de Gennes (BdG) Hamiltonian for the static DW, which has 
the following form: 

The self-energy function E describes the effects of the DW displacement. At the Gaussian level, we keep only the 
terms of the first and second order in X in its expansion: E = Ei + E2, where Ei = —XT-Li, with 

n.-(„\,^'t°) (10) 



V,AT 



and E2 ^ (XV2)-H2, with 

■Ho 



V^At 



^0 



Subsequent calculations are facilitated by two identities: 

mi\3) = -{E.^E,)W,,, (11) 
(*l^2b) = Y.^E, + E, - 2Ek)W,kWkj, (12) 

k 

where Wij = {i\'Vx\j) are the matrix elements of the generator of the DW translations. The first identity follows 
immediately from the expressions Hi = [Va;,Ho] and (j|[V2;, 'Hollj) = {Ej — Ei){i\Wx\j)- To prove the second identity, 
we observe that ■H2 = [V^;, [V2;,'Ho]]. On the other hand, 

(*|[V„[V,,Ho]]|j) 

= (E, + Ej){i\Vl\j} ~ 2(i|V,HoV,|j) 

= {E, + Ej) ^ W.kWkj - 2 EkW^kWkj. 
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The last line here is obtained using the completeness of the eigenfunctions. 

Inserting Eqs. ([5]), ([TT|). and ([T^ into the expansion of TrlnC/^^, it is straightforward to show that the terms linear 
in X vanish, while the quadratic terms can be collected into the expression (j4]) for the Gaussian action. The kernel 
is given by /C = /Ci + /C2, where 



/C2 = -Tr(S,5o) = -2^1^1^-^|M^. 



Since /C2 = — A^i(i^m = 0), the kernel vanishes for a stationary DW, as expected. Calculating the fermionic Matsubara 
sums, we arrive at the following result: 



,^1 ^ iE,-E,)[f{E,)-f{E,)] 



4 ^ -i^j 



where /(i?) = (e^^ + 1)^^ is the Fermi function. It is easy to see that the transitions between the BdG eigenstates 
corresponding to the same energy do not contribute to K.{vm). 

We note that one could also arrive at Eqs. (U) and (IT^ via a somewhat shorter route, using a change of coordinates 
X — X{t) — x' , y — y' , r — r', to transform into the reference frame co-moving with the DW. In this way, the 
invariance of the action under a static displacement of the DW is manifest from the beginning. A drawback of this 
approach is that the abovementioned transformation implies periodic boundary conditions for the fermionic wave 
functions as well as for the order parameter, which are actually inconsistent with a single DW. 

It is convenient to represent the Gaussian kernel (jl3p in the following form: 

l^(vra)^— / dede — 2 A/ e,e), (14 

where A/'(e, e') = Y.^J W^j\'^5{e - E^)S{e' - Ej), satisfying 7V(e,e') = N{e',e). To avoid dealing with ill-defined 
expressions for the matrix elements involving the bulk quasiparticle states, see Sec. llVl below. one can use the idenity 
(ITTt to obtain 

■^(^' = -uh^ E K*l^lb')l'^(^ - ~ E,), (15) 

ij 

where the summation is performed over the pairs of the eigenstates of "Ho, satisfying Ei 7^ Ej. The matrix elements 
of T-Li are well defined, because the order parameter derivatives are nonzero only in the vicinity of the DW. 



IV. QUASIPARTICLE SPECTRUM 

In order to calculate A/'(e, e'), we need to know the quasiparticle spectrum for the static DW. In the absence 
of external magnetic field, the 4x4 BdG Hamiltonian © can be written as a direct sum of two identical 2x2 
Hamiltonians, labelled by the spin projection a ='\, J,. From this point on we drop the spin index, restoring it only in 
the final expressions. For a DW parallel to the y-axis, the two-component wave function for each spin projection can 
be written as e*''''^4'(x), where '^{x) satisfies the equation 



^2 



2m* 
V AS(x) 



Ao(a;) 
2m* 



(16) 



Here fco 



kl and Ao(a;) = r]i{x){k^/kF) + V2{x){ky/kp). 



Since the superconducting order parameter varies slowly on the scale of the inverse Fermi wavevector, one can use the 
semiclassical, or Andreev, approximation"^^ and seek the quasiparticle wave functions in the form '^/{x) = e''^^^'0(a;). 
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where kx = ifeg. The slowly-varying envelope function ip = {u,v)'^ has the electron-like (u) and hole-like 
components, which are found by solving the Andreev equation 



= El}}. (17) 



The Fermi wavevector kp = (fcx, ky) — fci?(cos0, sin 6) defines the direction of semiclassical propagation of quasipar- 
ticles, along which the DW order parameter is given by 

Afcp(x) = 77i(x) cos6' -I- 772(2;) sin6', (18) 

and = vp cos9 [vp = kp/m*). Different models for the DW structure, see Sec. [TTl result in different semiclassical 
order parameters. However, the asymptotic values of A^^ (x) are fixed as follows: 

A_=Afc^(2;«-ed)=Aoe'^ , . 

A+ = Afc, [x » Cd) = Aoe'^e-^', ^'-"^ 

according to Eq. ([2]). 

Since A^^pix) — —Akp{x), the Andreev Hamiltonian defined by Eq. (IT71) has the property H^^p = —H^p. 
Therefore, the quasiparticle spectrum is electron-hole symmetric, in the following sense: if E is an eigenvalue of -ff^^, 
with the eigenfunction given by "0, then —E is an eigenvalue of H-k^, with the same eigenfunction. 

At given kp, the spectrum of Eq. ([T7|) consists of a continuum of scattering states in the bulk, with \E\ > Aq, and 
also discrete bound states (called the Andreev bound states, or ABS's) with the energies \E\ < Aq. Below we discuss 
general properties of the spectrum which are independent of a particular choice for the DW structure. 



A. Scattering states 



At each energy satisfying \E\ > Aq, there are two scattering states, labelled by L and R, corresponding to the two 
possible directions of propagation of the incident Andreev modes. Their wave functions can be found explicitly only 
far from the DW, i.e. at |a;| ^ ^d, where the order parameter is uniform. We have 



i^E,L{x) = C ■ 



4'E,B.{x) = C 



for the left-incident states and 

L+^e-^«"+ri^4"^^e^9^ a; » ^d 
for the right-incident states. Here q{E) = \J E"^ ~ t^j\vp^x\ > and 



(20) 



(21) 



(±) 



1 
71 



/ A± , vp^^q . \ 



V 



1 



VF,xq 

E 



I 



The Andreev scattering states are normalized, in the (S-function sense, with the normalization coefficient given by 

1 \ 



C{E) 



\E\ 



(22) 



The proof is presented in Appendix [B] 

The Andreev reflection and transmission coefficients j^(i?) and iL,_R(£') can be found by matching the asymptotics 
(|20|) and (|2ip and the solutions of Eq. p7|) near the DW. They satisfy the following general properties, independent 
on the details of the DW structure: 



I^LP + kip = ItflP + kflp = 1, t*p^rL + r^iL = 0, 
tL _ A- 
tR A+ ■ 



(23) 
(24) 
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It follows from Eqs. (l23l) and (|24l) that It^l = = t and Ir^l = |r/j| = r = ^/l — t^, and also that the scattering 
matrix defined as 

(25) 



is unitary. 

Expressions (1^^ and ([M)) follow from certain "conservation laws" for the Andreev equation. Let ipE,pi — 
{'UE.pi,VE,pi)'^ and '0£,p2 = {'^E,p2iVE,p2)^ be two solutions of Eq. ([T7| corresponding to the same energy, with 
Pi^2 — L 01 R. We define an analog of the Wronskian as follows: w[ipE,pi,i^E,p2] — pi<5'3V'-E,p2)- I* is straight- 

forward to show that dw/dx = 0, i.e. w['il}E,pn V'-E,p2] does not depend on x. Also, ip = aiip* corresponds to the same 
energy as ip and w[ipE,pi,ipE,p2] does not depend on x either. Therefore, 



'-^E,pi {x)uE^p2 (x) - ^B.pi (a;)wB,p2 (x) = const, (26) 
VE,pi {x)ue,p2 {x) - UE,pi {x)ve,p2 (x) = const. (27) 



The constants on the right-hand side can be calculated far from the DW, using the asymptotic expressions (I^Uj) and 
(CT . The properties and jM]) are obtained from Eqs. and ^7^, respectively. 

Explicit analytical expressions for the reflection and transmission coefficients can only be derived in some simple 
cases, in particular, for a "sharp DW" model, in which the DW thickness is sent to zero, see Appendix O These 
expressions can be used to find the asymptotics of t^^fi and r^^ji for an arbitrary DW of a finite thickness S,d at the 
energies close to the bulk gap edge, when the wavelength of the Andreev modes is much greater than ^d- It follows 
from Eq. jCU that 

tLM = 0, rLM = -l, (28) 

at \E\ = Aq. 

It is also possible to find the asymptotics of the reflection and transmission coefficients for an arbitrary DW at 
\E\ ^ Aq. At large energies, one can neglect the off-diagonal terms in the Andreev equations pT)). Then, the solutions 
for the left- and right-incident modes have the form 



^/2tt\vf,x\ V " / 



if E/vf,x > 0, and 



1 ( 

1 



^-iEx/vf 



W2t:\vf,x\ \ ^ / 



-iEx/vp_^ ^(^^ _ ( 1 piEx/vF^^ 



^y2^T\vF,x\ V / ' \/2n\vF,x 

if E/vf.x < 0. Comparing these expressions with Eqs. (PU)) and (HI]), we obtain: 

tL = tR ^ I, VL ^ rji = 0, for E/vF,x > 0, 

^+ 

A+ 

tL = 1, tii^ rL = TB.^ 0, for E/vF,x < 0, 

at \E\ > Aq. 



(29) 



B. Bound states 



At subgap energies, i.e. at \E\ < Aq, quasiparticles cannot propagate in the bulk, but Eq. (flT)) has solutions which 
are localized near the DW. It turns out that the number of such solutions can be expressed in terms of the properties 
of the scattering states: 

i , det 5(00) det 5(— oo) , , 

^- = ^+2^^" det5(Ao)det5Uo) ' ^''^ 

where S{E) is the scattering matrix defined by Eq. The proof is presented in Appendix [Dl The expression ([50)1 

plays the role of Levinson's theorem for the Andreev equation (recall that Levinson's theorem relates the number of 
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bound states of the Schrodinger equation to the phase shifts of the scattering states, see Ref. [3^. Note that there is 
a formal similarity between the Andreev Hamiltonian for a superconducting DW and the Dirac Hamiltonian in one 
dimension. The analogs of Levinson's theorem for the Dirac equation have been extensively studied, in particular, in 
the context of soliton charge fractionalization.— 

The determinant of the 5-matrix is a pure phase, which allows one to write the second term on the right-hand 
side of Eq. (|30l) in a more transparent form. From the asymptotics ([SS)) and ([29| . we have det5(±Ao) = — 1 and 
det S{±oo) = 6=*='^, where C = (26' — x) signu^.x- Therefore, ln[det S'(±Ao) / det 5(±oo)] = ztin =F iC + 2iTTn± and 

Nb ^l + n++n^, (31) 

where the integers n± are the winding numbers picked up by the phase of det S{E) as the energy varies between the 
bulk gap edges and the infinities. In particular, for the sharp DW model the determinant of the ^-matrix is given by 
Eq. (jCSp . from which we obtain = n_ — 0, therefore, Nb = 1, in agreement with the direct calculation of the 
bound states in Appendix [C] Different models for the DW structure yield different results for the ABS spectrum: 
while Nb = 1 for a sharp DW, one can have Nb > 1 for a DW of finite width. ^ 

At given kp = kp (cos 6, sin 6), the bound states have the energies Ea (a = 0,...,Nb — 1) and are asymptotically 

described by the wave functions ipaix) ~ g-V'^o-^i\^\/\^p,''\ ^ at \x\ 3> ^d- As the angle 6 varies between and 27r, 
the energies change, forming the ABS bands Ea{9). The BdG electron- hole symmetry is manifested in the following 
property: 

Ea{e)^~Ea{e + ^). (32) 

We use the index a = to label the branch of the ABS's whose energy vanishes at some 9 (zero modes). Note that if 
Eo{9) — then, according to Eq. (1321) . Eo{9 + Tr) =0 as well. The existence of zero modes is dictated by a topological 
argument, which relates the number of such modes with the difference between the fc-space topological invariants of 
the chiral order parameters in the two domains <^ Taking the spin into account, there are two pairs of spin-degenerate 
zero modes corresponding to the opposite directions of semiclassical propagation, i.e. four zero modes altogether. 
This is confirmed by the explicit calculation for a sharp DW model in Appendix [C] In general, fermion zero modes 
are present on any interface separating two superfluid or superconducting states with different topological charges, 
see Ref. [sS for a review. 

An important property of the ABS bands is the absence of degeneracies: 

Ea{9)^E,{9), (33) 

at any a and b. This can be shown as follows. Suppose that at some 9 there are two bound-state solutions of Eq. 
(ITTt . V'l and -02, corresponding to the same energy E. One can use Eq. (|27| with the constant on the right-hand side 
equal to zero (due to the exponential decay of the bound states at infinity) and obtain ipii^) = F{x)'4'2{x), where 
F(x) is a scalar function. Inserting this into Eq. (fT7|) . we have vp^x'^xF = 0, therefore F{x) = const, i.e. ipi and ip2 
in fact describe the same state. The ABS's can only become degenerate when vp^^ — 0, but at such directions of kp 
the Andreev approximation is not applicable. 

The bound states are responsible for a nonzero density of states at \E\ < Aq, which affects the system's low- 
temperature thermodynamics (for instance, there is a linear in T contribution to the specific heat, whose magnitude 
is proportional to the volume fraction occupied by the DWs, see Ref. [Til ), and also infiuence the Josephson current 
between two chiral superconductors.— 



V. DOMAIN WALL FRICTION AND MASS 



Now we are in the position to calculate Af{e,e') in the Gaussian kernel (fT4)) . In the semiclassical approximation, 
the states i and j in Eq. (|15|) correspond to the same Fermi wavevector kp, so that the matrix elements can be taken 
between the solutions of the Andreev equation (|17l) . We have 

+(p2(e,e')e(|e|-Ao)0(Ao-|e'|) 
+^2{e',e)9{Ao-\e\)9{\e'\-Ao) 
+^,{e,e')9{Ao-\e\)9{Ao-\e'\)]. (34) 

At given kp, Af{e,e') contains contributions from the transitions between different scattering states (ipi), between 
the scattering states and the bound states {(P2), and also between different bound states {(ps). The total intensity of 
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the quasiparticle transitions is obtained by summing over ah directions of semiclassical propagation and over the two 
spin projections: 



Here 



(pi(e,e') = 2^^|(e,p|g|e',p')P, 

kjp p.p' 

^2(e,0 = 2^^^ |(e,p|Q|6)p<5(e' - E,), 

kp P b 



(35) 



Q = 



V 





dx 



dAkj, jx) 
dx 




V 



with Akp(x) given by Eq. (llSp . p,p' = L or i?, and a,b = 0, ...,Nb — 1- The summation over kp amounts to the 
integration over the angle 0: 



p2n 



where Np = m*/27r is the Fermi-level density of states. It follows from the electron-hole symmetry of the quasiparticle 
spectrum that N{—e, — e') = A/'(e, e') and, therefore, V3i(— e, — e') = (^^^(e, e'). 

To calculate the matrix elements in Eq. (1351) one needs to know the wave functions of the bound and scattering 
states. Lacking such knowledge for a general DW, one can still make progress using dimensional arguments. Since 
the only energy and length scales of the problem are given by Aq and respectively, it is not difhcult to show that 



_A2 

V 



p b 



A2 



e t' 



(36) 



J2\{a\Q\b)\^5{e-Ea)Sie' -Ek 

ayib 



e e' 

a;;' a;;'' 



where Fi{x,x' ;9) are dimensionless functions. Inserting these expressions in Eq. psp . we obtain: 



2NfAI 



vf 



dO 



I cos 6*1 ' 
_ NfAI 

VF 



e e' 

Ao' Ao; ' 



(37) 



where fi are dimensionless functions. The Fermi surface angular integrals logarithmically diverge at — >■ ±7r/2 
(i.e. for the quasiparticles moving almost parallel to the DW) and have to be cut off at | cos 6*1 ~ a/AoTef ^ 1- 
In this way we obtain, with logarithmic accuracy: fi = ln(ei?/Ao)/i, where = kp/2m is the Fermi energy and 
Mx,x') = 2F,{x,x';TT/2) + 2F,{x,x';-n/2). Therefore, 



<P»(e,e') 



Vf 



Ao 



Ao' Ao 



(38) 



The expressions above can be calculated explicitly in the case of a sharp DWJ^ It follows from Eq. (jCip that 



Q = Ao 



p 

p* 



S{x) 
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where p = (A+ - A_)/Ao. Using Eqs. and (HI]), we obtain: 



■ sm 



X\ \xx'\t^{x)t^{x') 
2) V^^2—i^/^^r^=j 



1 



■ cos 



(39) 



Here t is the absolute value of the Andreev transmission coefficient, see Eq. (jC2[) . for which we have 



a:2 -cos2(0-x/2) 



Using Eqs. (gH]), (ED), and we obtain: 

2 . , y\ \x\t^{x) 



F2{x,x'-e) = -sin2(0-|) 



1 



1 



■cos(26l - x) 



vr 



■ sign (cos 9) sin(20 — x) 



Ao 



(40) 



where the bound state energy Eq is given by Eq. (IC6p . Finally, F3 = 0, since there is only one bound state at each 
6. One can see that Fi are nonsingular near the bulk gap edge, i.e. at |a;|, \x'\ — )■ 1, because of the vanishing of 
the Andreev transmission coefficients. This property actually holds for a general DW, since the wavelength of the 
Andreev states near the bulk gap edge diverges, making it possible to neglect the DW width at \E\ — Aq. 
Eqs. ([M]) and (PO]) yield the following expressions for the energy dependence in Eq. 



fi{x,x') 



■ cos 



2 X \xx'\\/x'^ — ly/x^' 



1 



2 [a;2 - sin^(x/2)][x^2 - sin^(x/2)] 



1 



1 . 2X 

xx' 2 



h{x,x') = -cos^ - 



X \x\\J x"^ — 1 



h{x,x') 



0. 



2 — i 



x 



— sin — ) + S (x' 



■ X 

sm — 
2 



(41) 

(42) 
(43) 



In some exceptional cases the logarithmic approximation, see Eq. (|38)) . might be insufficient. For instance, it follows 
from Eqs. (|41l) - (|15|) that the functions /,; all vanish for a sharp DW with x = tt. In such cases one should use the 
more general Eq. (|37|. 



A. Results 



Now we turn to the calculation of the effective dynamic action for the DW, see Eq. (|T4)) . At finite temperatures, 
the most singular contribution to the action at iy,n comes from e' close to e, which means that one can keep 
only the first term in the expression p4p . The contribution from the transitions between the bound states can be 
neglected, because the ABS bands are non-degenerate, see Eq. ([55)1 . We have 



where e — e — e' . The integration over e can be extended to infinity, since we are only interested in the low-frequency 
limit. Inserting here ipi given by Eq. (l38l) . we obtain that the leading frequency dependence of the kernel (fT3|) is 
non-analytic: lC{vrn) — 'n\^m\, with the viscous friction coefficient given by 

„r,.*Mi„(it).(^). («, 



VF \Ao/ V ^ 
where 

$(y -f / hix^x). 

It is not possible to calculate the integral and obtain an analytical expression for $(2/). Even in the simplest case of 
a sharp DW, it follows from Eq. (|1T]) that 

V „„-2 X f°° dx a;2 - 1 



cosh^(?/x/2) x"^ — sin^(x/2) ' 



10 



which still cannot be evaluated in a closed form. 

At low temperatures, all we need is the asymptotics of at j/ ^ 1, which is given by ^{y) ^ e^^. Therefore, 



at T ^ Aq. Physically, the DW friction is caused by the transitions between the bulk scattering states. These 
states absorb energy from the DW and then carry it away to dissipate into the thermal reservoir. Since the bulk 
quasiparticles are gapped, the temperature dependence of the friction coefhcient is exponential. 

Note that the expression (j44[) has the same general order of magnitude, but a different temperature dependence, 
as the friction coefficient of an ^ — i? interface in superfluid "^He (Refs. [20| and I21I). In the latter case, there are 
low-energy quasiparticles in the bulk with the momenta close to the A-phase gap nodes. The interface friction at low 
temperatures is dominated by the Andreev reflection of such quasiparticles off the interface, leading to a power-law 
behaviour r]{T) oc or T^, depending on the orientation of the orbital vector I relative to the interface. 

At r — 7> 0, the friction is negligibly small and the DW dynamics is dominated by inertia. Setting = inside the 
integral in Eq. p^ . we obtain /C(t'm) = where 

M^rdeTde'^ (45) 



Jo J -00 c — c 

is the effective mass per unit length of the DW. Inserting here Eqs. and we obtain 

M^C^^nfiJL), (46, 



VF VA 



^0 



where 



Ji J -00 [x-x'Y 7i 7_i [x-x'Y 7o [x-x'Y 

is a dimensionless coefficient. We see that the transitions between the electron and hole branches of the continuous 
spectrum (the first term on the right-hand side) , between the bound and the scattering states (the second term) , and 
between different bound states (the third term) all contribute to the DW inertial mass. 

Since Np = m* jlix in two dimensions, we have M ^ {m* /^o) ln(eF/Ao), where ~ w_f/Ao is the pair coherence 
length. At zero temperature ~ ^oj therefore, the ratio of the DW mass to the total fermionic mass per unit length 
contained in a strip of width is of the order of (Ap/ej?)^ ln(ef /An)^ i.e. very small. This is consistent with the 
estimate of the effective mass of the A — B interface: according to Ref. [20!, it is smaller by a factor of (Ao/ej?)^ than 
the total mass of the superfluid in the interface region. 



VI. CONCLUSIONS 



We have calculated the viscous friction coefficient and the zero-temperature effective mass of a domain wall in a 
chiral p-wave superconductor. The origin of both can be traced to the interactions of the DW with the Bogoliubov 
fermionic quasiparticles. The viscous friction is caused by the transitions among the bulk quasiparticle states induced 
by the DW motion. The friction coefficient is exponentially small at low temperatures, due to the bulk quasiparticles 
requiring thermal activation. The effective mass is determined by the transitions involving both the scattering states 
in the bulk and the Andreev bound states localized near the wall. 

The classical equation of motion for a DW can be written as MX + rjX — F, where the effective mass M and 
the friction coefficient r/ are given by Eqs. (|46p and pij) . respectively. The right-hand side contains the external 
driving force, which can come from the interaction of the orbital magnetization of the Cooper pairs with the external 
magnetic field, ■^'^ or from the coupling of the superconducting order parameter with the lattice deformation created 
by a sound wave<^ The latter mechanism offers a direct way of measuring the dynamical characteristics of the DW's 
by probing their contribution to the ultrasound attenuationr'^^ 

The DW dynamics for small deviations from equilibrium are determined by the Bogoliubov quasiparticle spectrum in 
the presence of a static DW. We presented a detailed investigation of this spectrum, including the general properties 
of the scattering states as well as the derivation of an analog of Levinson's theorem counting the number of the 
Andreev bound states. We did most of our calculations for a general DW with an arbitrary phase difference between 
the domains, without relying on any particular model for a DW structure. To illustrate the general formulas, we 
discussed in detail the case of a DW of zero width, for which one can make considerable analytical progress. 
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Our results are immediately applicable to clean neutral fermionic superfluids, such as '^He or cold atomic Fermi 
gases. In real superconductors, one has to take into account magnetic fields and screening currents, as well as disorder. 
In particular, impurities are expected to change the DW dynamics qualitatively, due to pinning. 
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Appendix A: Ginzburg-Landau description 

To develop a phenomenological understanding of the static DW structure, in particular, the origin of the phase 
difference between the domains, we use the Ginzburg-Landau (GL) free energy density given by _F = Fu + Fg, where 

Fu = a{\m\' + \V2r) + /3i(hiP + 1^21')' + + %T (Al) 

is the uniform part and 

Fg = Xi(V,77,)*(V,7y,) + K2iV^mnVjVj) + if3(V,?7jr (VjTyO (A2) 

is the gradient part. To avoid unnecessary complications we use the expression appropriate for the isotropic case. 
The chiral states q — Ao(l, ±i), with Aq = ^/\a\J4Pl correspond to the minimum of F„ if /3i, /32 > 0. 
The DW structure, see Eq. ([1]), can be written as follows: 

mix) ^ Aofi{x)e'^^^\ ry2(x) = Ao/2(a;)e'*(^)-^^(^), (A3) 

where /i^2 are dimensionless amplitudes of the order parameter components. The order parameter asymptotics are 
given by Eq. ([2]) . 

The origin of a nonzero phase difference x can be traced to the condition of vanishing supercurrent across the DW. 
The latter is obtained from Eq. (|A2p in the standard manner, with the result ji = 2lm{Kiri*\7ir]j + K^ri*^ jiy + 
K^rfjV jrii). Inserting here Eq. (jASp . we obtain: 

= 2A2(A'i23/f + ifi/|)(V,0) - 2XiA2/|(V,7) (A4) 

where A123 — Ki + A2 + K^. The presence in this expression of both the common and the relative phase gradients 
reflects the intimate coupling of the gauge and the internal (orbital) degrees of freedom in p-wave fermionic superfluids. 

Because of the current conservation we have V xjx — 0, therefore jx — const. The value of the transverse current 
is fixed by external sources and one can set jx = at all x. Then, Eq. (jA4[) yields a linear relation between the 
gradients of (j) and 7, which allows one to eliminate the common phase from the GL energy functional. The result is 
as follows: 

Fu = a^lifl + /I) + /3iA4(/2 + /2)2 + ^2A4(/4 + fl + 2/2/1 COS 27), 

T<C K f2 f2 (A5) 
Fg = ifi23A2(V./i)2 + K,Al{Vxf2? + lP7^Z%%2 Ag(V.7)'. 

Variational minimization of these expressions yields a system of three coupled nonlinear differential equations for 
fi,2{x) and 7(2;), subject to the boundary conditions /i^2(±oo) = 1 and 7(±oo) = ±7r/2. Using the solution of these 
equations, we can calculate the parameter x in Eq. ([2]): 

The value of x is manifestly non-universal, in the sense that it depends on the microscopic details. We note that, while 
the "locking" between (j) and 7, which results in the relation (jA6l) . is due to the condition jx = 0, the supercurrent 
along the DW remains nonzero. 

Due to the complexity of the equations for fi.2{x) and 7(0;), there is no exact analytical solution for the DW 
structure. To make progress, one can use, e.g. a constant-amplitude ansatz for the order parameter components;^ 
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which amounts to putting fi.2{x) = 1 at all x. Then, we obtain from Eq. (jA5[) the following expression for the free 
energy: 

F -(...) + i^A2(V,7)2 + 2/32 cos 2j, 

where the first term contains the 7-independent contributions and K — KiKi23/{Ki23+Ki). The variational equation 
for the relative phase has the form of a sine-Gordon equation, with a kink-like solution sin7(x) = tanh(a:/^d), where 



£,d = y^/4/32Aq is of the order of the GL correlation length and has the meaning of the DW thickness. From Eq. 
^61) we have 

In the weak coupling model, Ki — K2 — K3 (Ref. [H), therefore, x — ^/^- 

Appendix B: Normalization of the Andreev scattering states 

Let us consider two solutions, V'-E.p S'^.d '4'E',p', of Eq. (jl7p . corresponding to energies E and £", with p,p' — L 
or R. It is easy to show that ivF,x^ x tr{tp^^ ^astljE' ,p') = {E — £") tr(?/'^ ^V-'B'.pOj where "tr" denotes a 2 x 2 matrix 
trace in the electron-hole space (setting E — E' , we recover the "conservation law" for w['ipE,p, V'-E.p'Ij see Sec. IIV Ap . 
After integration, we arrive at the following useful identity: 

(E-E) tT{ipl,^j^'iljE',p')dx = ivF,x tT{tpl,^j^a3ipE',p') , (Bl) 

J Xl ' ' Xl 

which is valid for arbitrary Xi and X2 ■ 

The next step is to put xi = —i/2, X2 — and take the limit £ — > 00. The integral on the left-hand side of Eq. 
(|Bip becomes the inner product of the states iI^e,p and ipE'y, denoted by {E,p\E' ,p'). To prove the normalization, 
it is sufficient to consider the case E' — >■ E. The normalization integral for the scattering states of the Schrodinger 
equation contains a Dirac J-functionj^^ i.e. should be interpreted as a generalized function, and we expect the same 
to hold for the Andreev scattering states as well. After the substitution of the asymptotic expressions (pil)) and (|21l) . 
the right-hand side of Eq. (|B1[) contains the terms proportional to e^*''^"'"' ^^/^ and e^'*^'~'^ where q' — q{E'). 
The former terms oscillate fast at £ — )• 00 and can be neglected. 

Since the terms containing e^*^'^"'' ^^/^ oscillate fast unless q — q' , one can put E = E' in the pre-exponential 
coefficient. Then, using the properties (P^ . we obtain 

{E^ L\E'. L) . iE, R\E'. R) . ^^C%E) lim (B2) 
{E,L\E',R) = {E,R\E',L) = 0. 

The identity lim£_j.oo sin(a::^)/x = tt6{x) allows one to write the last factor on the right-hand side of Eq. (|B2p in the 
form 

^^^Mg^mi=^signES{E-E'). 
^— >oo E — E' 

Substituting in Eq. (|B2p the expressions for the normalization coefficients, we finally obtain: {E,p\E',p') = 
5pp,5{E-E'). 

Appendix C: "Sharp DW" model 

Many qualitative features of the DW quasiparticle spectrum can be illustrated using a simple model, in which there 
is a sharp boundary at a; = between the two domains with uniform order parameters of opposite chirality. The 
order parameter in the Andreev equation (I17p has the form 

Afc^(a;) = A_0(-a:)+A+0(x), (CI) 
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FIG. 1: The ABS energy for a sharp DW as a function of the direction of the quasiparticle propagation, for different values of 
the phase difference across the DW; X = (bottom panel), x — ""/^ (middle panel), and x = ^ (top panel). 



where A± are given by Eq. and d{x) is the Heaviside step function. The boundary condition for the Andreev 
wave functions is i^{+0) — ip{—0). 

For the scattering states {\E\ > Aq), one can use Eqs. ([20| and (|2T|) at all x. A straightforward calculation produces 
the following expressions for the transmission and reflection coefficients, 



tL = -: — tii = 



TL =rR = 



A+ (A+ + A_)wF,x9+(A+-A_)£;' 
Ao(A+ - A_) signi; 



(C2) 



(A+ + A^)vF,c,q + (A+ - A^)E' 
and also for the determinant of the 5- matrix, see Eq. ([25)1 : 



det5(i.) = ^^3±^, (C3) 



where A — tan(0 — x/2) signw^.a;. 

For the subgap bound states {\E\ < Aq), the normalized wave function has the form 



A. 



V^(^) = WTH^e-"'"'^'''"'^' I ET^^^signvF,, I, (C4) 



2\vF,x\ 



1 



where fl — ^y Aq ~ E'^. The upper (lower) sign corresponds to a: > {x < 0). Matching the wave functions at a; = 
we arrive at the characteristic equation 

E + zf^sign^F^^ _ A_ ^^^^ 
E — iflsignvF,x A+ 

It has only one solution, which can be found as follows. Let us introduce E = _E signfi?_a;. Since i?^ + fi^ = Ag, one 
can write E = Aq cos 6, n = Aq sin 9. It follows from Eq. jCS]) that e^*® = gK^e-x)^ therefore, 9 = 6'- x/2 + ttu 
and E = Ao(— 1)"" cos(6' — x/2). The integer n is found from the condition that > 0, which yields sign (sin 9) = 
(— 1)" sign [sin {6 — x/2)] = 1. Collecting all pieces together, we obtain the following expression for the ABS energy: 

Eo{e)=Aosie)cos(0-^), (C6) 
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where s{9) = sign [sin (0 — x/2) cos0]. 

Expressio n (IC6I) is valid for an arbitrary phase difference across the DW, thus generahzing the resuhs of Refs. 
andllSl. in which the bound states were studied for the (1, ±i) (i.e. x = 0) and (±1, i) (i.e. x = DWs. Note 
that the ABS energy is not a continuous function of 9, in generaL The discontinuities occur at the special directions 
of semiclassical propagation: a.t 9 — ±7r/2, i.e. for the quasiparticles moving parallel to the DW (in this case the 
Andreev approximation is actually not applicable and a more accurate treatment is needed), and also at 9 = x/2 and 
6* = x/2 + Ti", for which the DW is "invisible" to the quasiparticles, because A+ = A_. 

In Fig. [TJ we plotted the ABS energy for several values of x- In particular, if x = tt, then Eq. (IC6[) yields 
E = — Aosin^ = —Aoky/kp (see also Ref. [ll|): vanishing at ky = 0. The presence of zero modes is in fact generic: 
the ABS energy vanishes at = (x ± 7r)/2, resulting in low-energy quasiparticles bound to the DW. Taking the spin 
into account, we have two pairs of spin-degenerate zero mode branches. 

We would like to note that one can also use Eq. (|C5[) to obtain the ABS spectrum for a DW between two isotropic 
s-wave superconductors, which can be realized as a Josephson junction with the phase difference x- In this case, the 
right-hand side of Eq. ()C5| is equal to e~^^ and the bound state energy is given hy E = —Aq signwF.a; cos(x/2) (Ref. 



Appendix D: Proof of Eq. (fSO]) 

We start with the Andreev Hamiltonian H defined by Eq. P7)l . with the order parameter written in the form 
= |Afcy|e''^, where the amplitude and the phase have the following asymptotics: ip{x) = 6* at x <C —£,d, 
fi^) = X ~ ^ at a; ^d, and |Afcp(a;)| = Aq at |a;| ^d- To represent the DW as a locahzed perturbation, 
we remove the phase from the off-diagonal elements by a unitary transformation as follows: U^HU = H, where 
U = e"^'^3/^. The transformed Hamiltonian is given hy H = Hq + SH, where Hq = —ivp.x^a^ x + AqCTi describes the 
Bogoliubov excitations in the uniform p-wave superconducting state, and 5H = vf,x'^' {x)(Tq/2 + (\Akp{x)\ — Ao)a'i 
is a perturbation which is nonzero only near the DW. 

The number of the bound states for H is the same as for H. To calculate the latter we observe that there is a 
one-to-one correspondence between the eigenstates of H and iJoj which are found from the equations Hipi = Eiipi and 
HQipf'' — El^^ip'f \ respectively (to verify this one can introduce H\ = Hq + X6H and consider the smooth evolution 
of the spectrum as the parameter A varies from to 1). The total number of states is "conserved" , which is formally 
expressed by the formula de[p{e) — po(e)] = 0, where 



p(e)=lim/ dxYS{e-E,)tTi)l{x)^,{x) (Dl) 

is the density of states for H and 



po(e)=lim/ dxY 8{e-Ef)\Yi,f^\x)^f\x) (D2) 



is the density of states for Hq. One can write p = pb+ PSi where pB [ps) is the contribution of the bound (scattering) 
states. On the other hand, po is nonzero only at \E\ > Aq. Therefore, 

= / pBde+ (ps - Po)de ^ Nb + [ps ~ Po)de. 

J\E\<Aa J|£;|>Ao J\E\>Ao 



Inserting here Eqs. (|Dip and (|D2p and using the energy E and the direction of propagation p = L,R to label the 
scattering states, we arrive at the following expression for the number of the bound states: 



f di? hm ^ [trV;kp(^)^i^,p(^)- tr<'pt(x)<)p(a:)j (D3) 



p=L,R 



Since the eigenstates of H and H are related by a unitary transformation, ipE.p{x) = e "^'•'^■*'^^/^V'-B,p(a;), one can 
replace 'ipE,p{x) in Eq. (|D3p by ipE,p{x). 

The right-hand side of Eq. (ID3p is ill-defined because it contains the difference between two infinite (5-function) 
normalization integrals for the scattering states with and without the DW. To make sense of this expression, we split 

the energies of the eigenfunctions as follows: ipE^pipE,p - V'^'^Vb,], = ['^E,p'^E',p - '0s,pVb>]-B'^b- Then, one can 
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use the identity (jBl[) to express the integral over coordinates in Eq. (jD3|) in terms of the asymptotic values of the 
eigenfunctions far from the DW. In this way, we obtain 



Nb = -iVF,x 



dE lim P, 



(D4) 



Bl>Ao 



where 



P= lim 

E'^E E - E- 



p=L,R 



/2 ■ 

-1/2 



The asymptotics of '4'e,p are given by Eqs. (|20p and (|21l) . while the expressions for V'^p can be obtained from Eqs. 



20)) and (|2T|) by setting t^.i? = 1 and r^.n = 0. After some straightforward algebra, we have 

1 



lim 



E'^E E - E' 



A{E,E')e~''^''~'^'^^/^ + B{E,E')e''^'^+'i">^^^ + B*{E',E)e''^'^+''' 



)£/2 



(D5) 



where 



AiE,E') = [tl{E)tL{E')+rl{E)rR{E') - l]a++{E,E') - [rl{E)rL{E') + tl{E)tR{E') - l]a__(i?, £;'), 
B{E,E') = rR{E')a.+ {E,E') - rL{E')a+_{E,E'), 



and 



a3,.(S,£;') -C{E)C{E') 



1 + s 



, , ,VF.xq'\ . . w //t VF.xq\/^ ,VF,xq' 

1 + s'-^j Sign Sign £; - ^ (^1 - s-^ j (^1 - s'-^ 



with s, s' = ±. 

It is easy to see that A{E,E) = 0, due to Eq. and B{E,E) = 0. Therefore, 



P 



dA{E, E') 



dE' 



dB{E,E') 



E=E' 



dE' 



Jql 



dB*{E',E) 



E=E' 



dE' 



-iql 



(D6) 



E=E' 



In the first term, one has only to differentiate the coefficients in front of a++ and a , because |<i| = jt/fl = t and 

Vl\ = \i'r\ = VI — t^. In the last two terms, one has only to differentiate a |- and The result looks as follows: 

P ^ P^+ P2, where 



Pi=- 



1 



dtj 



drj 



dtR 



^-—signE^t,^+rl^+tl — 



' dE 



(D7) 



and 



P9 



^C^^Im 

VF,x \E\ 



In Pi, one can use the definition of the S'-matrix, Eq. p5|). to represent the expression in the brackets as tT{S^dS/dE). 
In P2, we observe that in the limit £ ^ 00 only small q are important, corresponding to the energies close to the buk 
gap edge. For these energies, the asymptotics (1^ hold and we have + rp ~2. Therefore, 



liniP.^^C^^hm^ 

1^00 vf.x \E\ 1^00 q 



VF,, 



+ Ao) + - Ao 



Inserting Eqs. (|D7l) and ((D8)) into Eq. (|D4|) . we obtain 



2ti 



dE sign tr ( S"^ 



B|>Ao 



.dS_ 
'dE 



(D8) 



(D9) 



The unity on the right-hand side originates from the states located exactly at the bulk gap edge, which give rise to 
the i5-functions in Eq. ()D8|) . When integrated over energy, each of these (5- functions contributes 1/2 to Nb- Finally, 
using tY{S'^dS/dE) = d{\nAetS)ldE, we arrive at Eq. 
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